Tuning odd triplet superconductivity by spin pumping 
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We study proximity effect in diffusive ferromagnet/normal- metal/superconductor junction with 
precessing magnetization of the ferromagnet. We find that the odd-frequency pairing induced in 
the normal metal is modified by spin pumping from the ferromagnet and hence can be tuned by 
changing the precessional frequency. At the frequency corresponding to twice the superconducting 
gap, the odd-frequency pairing is strongly enhanced. We find a crossover from the even- to the 
odd-frequency superconductivity in the normal metal by tuning the precessional frequency. This 
gives a clearcut signature of the odd-frequency superconductivity observable by scanning tunneling 
microscopy. According to the pairing symmetries in the normal metal, we find a crossover from the 
gap to the peak structure in the tunneling conductance between the normal metal and a scanning 
tunneling microscope tip. 



I. INTRODUCTION 

Generally, superconducting correlations can be even or 
odd in frequency, depending on their symmetry with re- 
spect to the time axis. In accordance with the fermionic 
statistics, even-frequency superconductors are character- 
ized by the spin-singlet even-parity or spin-triplet odd- 
parity pairing states, while odd- frequency superconduc- 
tors are grouped into the spin-singlet odd-parity or spin- 
triplet even-parity pairing states. 

The possibility of the odd-frequency pairing state in 
a uniform system was discussed in the literature^^^ 
although its realization in bulk materials is still 
controversial. The odd-frequency pairing state has 
recently been predicted in inhomogeneous super- 
conducting systems i&^&^ i 10 ' 11 ! 12 ' 13 In diffusive ferro- 
maget/superconductor junctions, odd- frequency pairings 
emerge due to the broken symmetry in spin space 
In the diffusive ferromagnet, only s-wave symmetry of 
the pair amplitude is allowed due to impurity scatter- 
ing and hence triplet pairing inevitably belongs to the 
odd- frequency superconductivity, which is also called odd 
triplet superconductivity. 

Recently, the interplay between ferromagnetism 
and superconductivity in ferromagnet/normal- 
metal/superconductor (F/N/S) junctions has been 
studied using spin-active boundary condition^ Another 
recent progress is the study of the Josephson effect in 
S/F/S junctions, where it has been established that the 
ferromagnetic spin dynamics play a crucial rolej 15 ' 16 ! 17 
The F/N/S proximity effect in the presence of spin 
dynamics, however, still remains to be fully understood. 
Moreover, although the long-range proximity effect by 
the generation of the odd-frequency triplet pairing has 
been observed in recent experiments^^ no experiment 
has succeeded in controlling the magnitude of the 
odd-frequency pairing to this date. 

In this paper, we propose an experimental setup which 
makes it possible to tune the magnitude of the odd- 
frequency pairing. We study the proximity effect in dif- 
fusive F/N/S junctions with precessing magnetization of 



the F layer. We find that the odd-frequency pairing in- 
duced in the N layer is modified by spin pumping from 
the F and therefore can be tuned by changing the preces- 
sional frequency. At the resonance frequency (i.e., twice 
the superconducting gap), the odd-frequency pairing is 
strongly enhanced. We find a crossover from the even- to 
odd- frequency superconductivity in the N layer by chang- 
ing the precessional frequency. This is reflected in the 
tunneling conductance between the normal metal and the 
scanning tunneling microscope (STM) tip as a crossover 
from the gap to the peak structure. 



II. FORMULATION 
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FIG. 1: (Color online) Schematic of the model for a 
ferromagnet/normal- metal/superconductor (F/N/S) junction 
with precessing magnetization of the F layer. A fictitious ex- 
change field is dynamically induced in the N layer, which in- 
terplays with the superconductivity induced by the proximity 
effect. The electronic structure in the N layer can be studied 
by scanning tunneling microscopy (STM). 

We consider a diffusive F/N/S junction with precessing 
magnetization of the F as shown in Fig. Q] The F/N 
interface is located at x — while the N/S interface is at 
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x = L. The spin relaxation in the junction is assumed 
to be weak. Then, a spin density is pumped into the N 
from F,— while the superconductivity is induced in the 
N by the proximity effect. Therefore, we can study the 
interplay between ferromagnetism and superconductivity 
in the N layer ^ 

Let us take the time-dependent exchange field in the 
F to be directed along 



m(i) = (sin 9 cos fii, sin 9 sin fit, cos ( 



(1) 



Here, Q is the precessional frequency around z axis and 
9 is a constant tilt angle. In the rotating frame, such 
precession can be viewed as a difference between spin- 
resolved chemical potentials along the z axis, and the 
noncollinear effective exchange fields in N and F can gen- 
erate a long-range proximity effect i 16 ' 21 ! 22 To investigate 
the proximity effect, we use the unitary transformation: 



g(t,t')^U\t)g(t,t')U(t'), 



(2) 



where U(t) = cxp (— transforms from the lab- 
oratory into the spin-rotating frame. Then, the prob- 
lem reduces to the stationary one^ In fact, we have 
W{t)u\{t) ■ crU{t) = m(0) • a- = m • a. Here, g is the 
retarded component of the quassiclassical Green's func- 
tion, which is 4x4 matrix in spin<g>Nambu space. We 
parameterize g as 



9 = t 3 ® (g + g • <r) + n <8> (f s +f t -<r), 



(3) 



where g and f s are scalars, and g and ft are three- 
dimensional vectors.— Oi and Tj (i = 0,1,2,3) are the 
unit and Pauli matrices in the spin and Nambu spaces, 
respectively, er = [ax, a-i, 03) is the vector of Pauli ma- 
trices. 

The Usadel equation in the N in this rotating frame 
has the formi^ 2 ^ 



DV(gVg) + [ist 3 + i(Cl/2)<r 3 ®T 3 ,g}=0, 



(4) 



where D is the diffusion constant and e is the quasi- 
particle energy measured from the Fermi level. See Ap- 
pendix [5] for derivation. To take into account the mag- 
netic proximity effect, we consider a low-transparency 
spin-active interface at the F/N contact, with the bound- 
ary condition on the N side given b y 14 > 26 



2jB^gd x g = [r 3 + r^fm ■ a) ® r 3 , g] 



(5) 



Here, jb = RbL/Rd£, and 7^ = G^/Gt, in terms of the 
interfacial resistance parameter Rf>, the diffusive resis- 
tance of the N Rd, the N coherence length £, the imag- 
inary part of the mixing conductance 2 ^ G^ and the in- 
terfacial tunneling conductance Gt ~ The mag- 
netic proximity effect is governed by 7^.— We will for 
simplicity disregard the real part of the mixing conduc- 
tance and the spin dependence of the F/N conductance, 
which should not affect our findings qualitatively. We 



make use of the well-known Kupriyanov-Lukichev bound- 
ary conditions 2 - for the N/S interface: 



^iBigd x g = \gs,g] , 



(6) 



with the bulk Green's functions gs- For simplicity, we 
assume the same jb parameter at both interfaces. 

Focusing on the tunneling regime with weak supercon- 
ducting correlations in N, we linearize the Usadel equa- 
tion with respect to the anomalous Green's function in 
the N. The linearized Usadel equation reads 



Dd 2 J s + 2ief s - 
Dd 2 J t + 2ief« 



2if t ■ h 
- 2if s h 



0, 
0, 



(7) 



with h = — Slz/2 and z = (0,0, 1). The general solution 
in the N is thus 
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with k± = ^2i(-e±Q./2)/D and k a = ^/~2ie/D. The 
coefficients A, A' , B, B', C and C are determined by the 
boundary conditions. Here, fx is the transverse triplet 
component, which is long ranged at low energies, and 
/ 3 is the longitudinal triplet component, which is short 
ranged, similarly to the singlet component f s & It follows 
from the boundary conditions that the second component 
of ft is zero. 

Linearizing the boundary conditions at the N/F inter- 
face, we have 



lBid x f s = f s + H<t>ft ■ m , 
iB^dJt = ft + «70/ s m ■ 



(9) 



These boundary conditions show that the presence of f s 
generates the triplet components with f t || m as long as 
70 =/= 0. Similarly, the linearized boundary conditions at 
the N/S interface have the form 



lBtd x f s + g° s f s + g 3 s f 3 = f° s , 
^d x f t +g s f t + g y s z = f s z. 



(10) 



Here, the bulk Green's functions in the S, gs and 
fs, are given by g s = (g s ^0 + 9s a 3) ® r 3 and f s = 
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(/s°"o + ft^s) ® T ij where (in a convenient gauge) 
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(11) 



After the unilluminating manipulations to solve Eqs. ©- 
([TO]) , we obtain the desired solution of the Usadel equa- 
tion. The explicit form of the solution is rather compli- 
cated and is given in Appendix IB1 

In general, the tilt angle 9 is determined by solving the 
Landau-Lifshitz-Gilbert equation: 



dm(t) 
dt 



-jm(t) x H e ff(t) + am(t) x 



dm(t) 
dt 



(12) 



where 7 is the gyromagnetic ratio, H e g- is the effective 
magnetic field and a is the Gilbert damping constant. 
Driving the ferromagnet by an applied transverse ro- 
tating magnetic field h v {, in the presence of an effec- 
tive static dc field -f/dc, we find a small transverse mag- 
netic field in the rotating frame of reference. This is 
disregarded in our analysis, assuming a <C 1, so that 
9 ~ hrf/aH^c 3> htf/Hdc at the ferromagnetic resonance, 
which is usually the case in realistic ferromagnets. We 
note that the spin-pumping effect is enhanced close to the 
ferromagnetic resonance, which should in practice be the 
optimal driving regime for the magnetic dynamics. Note 
that the tilt angle of the rotating magnetization 9 may in 
general be time dependent. In this paper, we focus on the 
regime when d\n9/dt <C (n, 9 2 /ri), so that the appropri- 
ate nonequilibrium spin state is fully developed before 
9 changes appreciably. Here, 1/t, is an effective spin- 
injection rate proportional to the mixing conductance 
at the F/N interfaced In addition, the spin-relaxation 



rate is assumed to satisfy l/r 8 f <C n, so that the spin 
memory is preserved during a cycle of precession and 
l/fsf 9 2 /ri, so that the developed spin accumulation 
does not decay. We will, therefore, fix 9 and disregard 
spin relaxation in the junction in the following. 

Next, let us consider tunneling current between the 
N and the STM tip. See Fig. [TJ Below, we consider 
two cases: spin relaxation in the STM tip is weak or 
strong, on the scale of its characteristic spin-injection 
rate. Weak (strong) spin relaxation effectively leads to 
a state of equilibrium in the STM tip, in the rotating 
(lab) frame of reference. To realize the situation that the 
state of equilibrium is effectively reached in the rotating 
frame of reference, the STM tip should be spin-polarized 
by the fictitious field, which would in practice require 
either connecting the precessing ferromagnet directly to 
the tip or using magnetized tip which itself is precessing 
in synchronization. Experimentally, this may be more 
difficult to realize than the strong spin-relaxation regime 
in the STM tip. (Note, however, that a spin-resolved 
STM technique has been recently developed^) One may 
furthermore think that the limit when n ~ A is problem- 
atic for the superconductor, since the gap should then be 
strongly suppressed. However, assuming some spin relax- 
ation in S, the out-of-equilibrium pumped spins do not 
cause significant Cooper-pair depairing in the supercon- 
ductor, while the effective spin relaxation in N is weak, 
so that the pumped spin distribution is unaffected in N. 

Now, let us first consider the case of the STM tip with 
weak spin relaxation. We define the 8x8 matrix Green's 
functions in STM and N, §stm and §n, as 



9STM 



T 3 9STM 







-T3 



9N 



9§ 9§ 
9n 







(13) 



Here, R, K and A represent the retarded, Keldysh and 
advanced components, respectively. The current between 
the N and the STM tip is given by 



cfeTr (r 3 [gsTM,9N 



K 



(14) 



where 



[gsTM,gN] K = 9stm9n - 9n9stm + t 39n + gN T 3 ■ (15) 
Let us apply a bias voltage V to the STM tip, so that we have 

9STM — T 3 (fsTM + fsTM T 3) + (fsTM + fsTM T 3) T 3 = %fsTM + ^fsTM T 3 i 

f s % Ia = [tanh{(e + el/ + crn/2)/2r} + (-)tanh{(e-ey-crn/2)/2r}]/2, (16) 
with a = ± for spins up and down, and = tanh(e/2T) (</^ — gfj). Therefore, we obtain 

Tr (r 3 [gsTM, 9n] K ) = £ 2/J TM , CT « CT - < CT ) = - £ 4/| TA/ ^ Re , (17) 
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and, hence, at zero temperature: 

dl/dV ~ l N ° (~ eV ~ an / 2 ) + N " ( eV + CT °/ 2 )] 

I 



(18) 



Here, N a = Reg^ denotes the density of states for spin 
a in the N. To calculate the spin-resolved density of states 
N a , we can use the relation g = —f s ft, which is obtained 
from the normalization condition g 2 = 1 applied to the 
parametrization ([3]). 

In the case of the STM tip with weak spin relaxation, 
the state in the STM is equilibrated in the rotating frame 
of reference, and we obtain the relevant results by setting 
fl = in Eqs. (fT6|) and (fl"8]) . The tunneling conductance 
then coincides exactly with the density of states. 



III. RESULTS 

In the numerical calculations, we introduce a small 
imaginary part of the energy, e — > s + iS, to regularize 
singularities. Physically, 5 captures an effective depairing 
rate for Cooper pairs. We choose the following param- 
eters: 7s = 100, 8 = tt/12, L/i = 10, A/E Th = 0.1, 
and S/E Th = 0.01. Here, E Th = D/L 2 is the Thouless 
energy. Typically, for L/£ = 10, E Th ~ 0.01 meV. 



8 =/= 0,7r/2 in this limit. Thus, we can expect that by 
tuning fl close to 2A one can control the magnitude of 
the odd-frequency superconductivity for sufficiently large 
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A. Weak spin relaxation in the STM tip 

We start by considering the case of weak spin relax- 
ation in the STM tip. We first discuss some limiting cases 
and give qualitative picture of our main results. Below, 
we focus on the low-energy limit, i.e., e — > 0, where the 
superconductivity is manifested most strongly. At e = 0, 
the bulk Green's functions in the S are given by 



9s = 0, g 3 s 
fl 



■in/2 



Va 2 - WW 



- ("/2) : 



, /l = o. 



for n < 2A and 



9°S 



n/2 ,,1 = 0, 



V(fi/2)2 - A 2 



ft 



0, fl 



iA 



- A 2 ' 



(19) 
(20) 

(21) 
(22) 



for n > 2A. In the limit of n — > 2A, these Green's 
functions diverge and the proximity effect is resonantly 
enhanced as seen from the boundary conditions at the 
N/S interfaces. Taking the limit 7^ — > 00 for the bound- 
ary condition at the F/N interface, we have f s — > 
and fi sin 8 + f 3 cos 8 —> 0. The singlet component van- 
ishes while the triplet components can remain finite for 



FIG. 2: (Color online) Anomalous Green's functions as a func- 
tion of il/ETh at the N/S interface [x = L) with (a) 7^ = 100 
and (b) 7^ = 10 and at the F/N interface (x = 0) with (c) 
70 = 100 and (d) 74, = 10. Here, we set e = 0. /1 and fs 
are the odd-frequency and / s is the even-frequency pairing 
amplitudes. 

Next, let us plot the anomalous Green's functions us- 
ing Eq. ([5]). Figure [5] shows anomalous Green's functions 
at the N/S interface (x — L) as a function of O/i^Th with 
(a) 70 = 100 and (b) 7^ = 10. Note that at e = 0, f\ and 
/3 are purely imaginary while / s is a real number. We 
find that a resonant peak appears at n = 2 A. As seen in 
Fig-IHa), for n < 2A, the singlet component dominates, 
while, for n > 2A, the triplet component dominates. 
This can be understood by the fact that for fl — > 2A — + , 
fs diverges while for fl — » 2 A + + , /J diverges. Thus, 
one can control a crossover from the even- to the odd- 
frequency superconductivity by changing fl, which is tun- 
able by the external magnetic field. When 7^ is reduced, 
the singlet component is enhanced as shown in Fig. [2b). 
We also show the anomalous Green's functions at the 
F/N interface (x — 0) as a function of fl/E^ in Fig.f^c) 
for 7^ = 100 and Fig. [^d) for 70 = 10. A peak also ap- 
pears at fl = 2A. Compared to the results at x = L, the 
long-range triplet component f\ has a large magnitude, 
which is controllable by tuning the frequency fl. 

To date, a hallmark of the odd-frequency pairing has 
been considered to be the long-ranged proximity effect in 
the presence of magnetism. 6 However, another aspect of 
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FIG. 3: (color online) The density of states normalized by its 
normal state value as a function of e/Erh with 7^ = 100 at 
(a) x = L and (b) x = 0. 



FIG. 4: (color online) The normalized tunneling conductance 
as a function of eV/ETh with 7^ = 100 at (a) x = L and (b) 
x = 0. 



this pairing has been recently appreciated: The density 
of states in the presence of the odd-frequency pairing is 
enhanced, acquiring a zero-energy peak within the gap 
struct ure£ i 21 ' 31 i 32 Using general relations for the conju- 
gate Green's functions^ 3 - we have that / s (e) = /*(— e) 
and ff(e) = — f t *(— e). Hence, we easily obtain g 2 
I — |/ s (0)| + I ft (0)| at e = 0, from the standard nor- 
malization condition g = 1. Therefore, the density of 
states, which is given by Reg, is enhanced by the gen- 
eration of odd- frequency pairing (ft) and suppressed by 
the presence of the even-frequency pairing (f s ) at e = 0. 

The density of states normalized by its normal state 
value, N — Re^/I — (f 2 + j\ + /f ), as a function of 
e/Erh is shown in Fig. [3] setting 7^ = 100 at (a) x = L 
and (b) x = 0. We find a crossover from the gap to the 
peak structure in the density of states with increasing 
f2/-&rh- This reflects a crossover from the even- to the 
odd- frequency superconductivity in the N, upon changing 
the precessional frequency. Note that, as shown above, 
the density of states exactly coincides with the tunnel- 
ing conductance between the N and the STM tip, in the 
weak spin-relaxation regime of the STM tip. 



B. Strong spin relaxation in the STM tip 

Here, we discuss the more realistic case with strong 
spin relaxation in the STM tip, so that the tip is equi- 
librated in the lab frame. Let us define ctt as dl/dV 
normalized by its normal state value. Figure U shows ot 
as a function of eV/ ETh for 7^ = 100 at (a) x — L and (b) 



x = 0. We see the gap structure at fi = 0. For larger Vt, 
the zero-bias conductance grows, and a sharp zero-bias 
peak emerges at f2/i?Th = 0.1. This can be understood 
as follows: Let us consider the term N^(eV + f2/2) in 
Eq. (|18p . for instance. Then, the argument e — 0/2 ap- 
pearing in Eqs. (j4|) and (fTT|) is effectively replaced by eV 
while e + fl/2 is replaced by eV + il. Similar replace- 
ments should be made for other terms in Eq. (|18p . Thus, 
coherence peaks will appear at eV = fi ± A. Therefore, 
when = A, a resonant peak appears at the zero bias. 
This can be also attributed to the emergence of the odd- 
frequency pairing, as shown in Sec. IIII Al Therefore, we 
find that even in the case of the strong spin relaxation 
in the STM, we can also find the evidence of the odd 
triplet superconductivity (zero-bias peak) tuned by the 
precession frequency. 



IV. CONCLUSIONS 

We have studied the proximity effect in diffusive 
F/N/S junctions with precessing magnetization of the 
F layer. We found that a tunable odd-frequency pairing 
in the N is governed by spin correlations in the N that 
are induced by the ferromagnetic precession. At the res- 
onance frequency, Q — * 2 A, the odd- frequency pairing 
is strongly enhanced. We unveiled a crossover from the 
even- to the odd-frequency superconductivity in the N 
by tuning the processional frequency. This is manifested 
as the crossover from the gap to the peak structure in 
the tunneling conductance between the N and the STM 
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tip, where we investigated two regimes of the weak and 
strong spin relaxation in the STM tip. 

Our results can be experimentally accessible and pro- 
vide a novel way to amplify the odd triplet superconduc- 
tivity using the magnetization precession in the F/N/S 
junctions. These characteristics should be observable by 
a scanning tunneling microscope or other tunneling ex- 
periments and, hence, may serve as a smoking gun to 
reveal the odd-frequency superconductivity. 

Recently, the F/N/S trilayer junctions have been fab- 
ricated to study behavior of the superconducting phase 
transition. 34 ' 35 Also, a ferromagnetic resonance experi- 
ment has been performed in an Nb/permalloy proximity 
system.- 3 ^ In the light of these advances, it appears re- 
alistic to verify our predictions experimentally by using, 
e.g., permalloy/Cu/Nb junctions. The predicted reso- 
nance frequency corresponds to the ^MHz range and is 
easily achievable by the present-day experimental tech- 



nique. It should be remarked, however, that the required 
condition for the spin-flip relaxation rate, l/r s f -C fi, is 
rather stringent in this low-frequency limit, requiring low 
temperatures and very clean normal interlayer. 



Acknowledgments 

We are grateful to A. Brataas, Ya. V. Fominov, 
M. Houzet, J. Linder and J. Inoue for helpful discussions. 
This work was supported by the JSPS (T.Y.) and by the 
Alfred P. Sloan Foundation (Y.T.). 



APPENDIX A: DERIVATION OF THE USADEL 
EQUATION 



J 



Here, we present the derivation of the Usadel equation in the presence of the exchange field. The matrix Green's 
function in particle-hole <E> spin space can be defined as£ 



G = -i 



T c a i a \ 
T c a \ a \ 



'T c a r a\ 

' T c a A 
\ T c a \ a \ 
\ T c a Wl] 



Tca\a, 



'T c a T a T ) \ 

T c a[a T 
Tca^aiJ j 



(Al) 



where Tq is the time-ordering operator along the Kcldysh time contour. This basis is obtained from the conventional 
basis, 



G 



by the following transformation: 



TcaA 
' Teat at 



Tea , a 



ri, 



'T c aja[ 
[ T c a Al, 



T c a|a T 
,T c a{a T 



G — > UGU\ 



Tca i a i 
T c a|a^ 
T c a\a^ J 



(A2) 



(A3) 



with 



U = 



1 

a x 



(A4) 



With the quasiclassical approximation, we obtain the quasiclassical Green's function g from G in this basis. 6 
Then, the pair potential is transformed as 



A 



Acr 2 
A*a 2 

A(7 2 C7l 

A*a x a 2 



(nRe A - r 2 Im A) ® a 2 



-iAa 3 
A* era 



(r 2 Re A + nlm A) <g) a 3 , 



(A5) 
(A6) 



and also we have 



a 

a* 



cr 

ox<y*o\ 



(to ® ax, T Q <8> cr 2 , r 3 <8> a 3 ) = S . 



(A7) 
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The Usadel equation is, therefore, transformed by U as 6 

DV(gVg) + [ist 3 - i(h ■ a) - (tiRc A - r 2 ImA) ® cr 2 ,g) = 
-> DV{gVg) + [ier 3 - i(h- S) - (Y 2 Re A + nImA) ® a 3 ,g] = 0. 

Next, we consider the following unitary transformation)^ 





" 7T 






7T 


exp 


l-T 3 6 


5 <J 3 


exp 


-i~a 3 



such that 

WiV 4 = (7i (8) t 3 , Fo- 2 F f = (7 2 O r 3 , Ut 2 <8> e^U 1 = ri , Vt x <8> cr 3 U f = -r 2 : 
employing the relation 

e ien<T = cos 6» + i sin 6>n ■ er , 

for an arbitrary unit vector n. The Usadel equation is finally transformed by V as^ 

DV(gVg) + [ieT 3 -i(h-S)-(T 2 ReA + TiImA)<g>(T 3 ,g] = 
-> DV{gVg) + [ier 3 - i(h ■ tr) ® t 3 - (r x Re A - r 2 Im A) <g> C7 , 5] = . 



(A8) 
(A9) 



(A10) 

(All) 
(A12) 



(A13) 
(A14) 



This representation of the Usadel equation for the transformed Green's function has the convenience of the explicit 
symmetry with respect to rotations of the exchange field h. 

Assuming that A is real, the equation for the / component yields 



DV(gVf + fVg) + 2ief - i(fh ■ tr + h • trf) = Ag - Ag . 
Linearizing this equation with respect to superconducting correlations, we find 

DV 2 f + 2ief - i(fh ■ a + h ■ trf) = -2A . 

Setting / = f s + ff ■ <7|2£ we have 

DV 2 (f s + f t ■ a) + 2ie(f s +f t -tr)- i(2f s h ■ tr + 2f t • h) 



2A. 



giving finally: 



DV 2 f s + 2ief s - 2if t h = -2A, 
DV 2 i t + 2ief t - 2if s h = 0. 



(A15) 
(A16) 
(A17) 



(A18) 
(A19) 



Here, ft and h are three-dimensional vectors. As shown in Ref. 28, f s and f t || h are short ranged while ft_Lh is long 
ranged, at low energies. 



APPENDIX B: SOLUTION OF THE USADEL 
EQUATION 



r 



A = -(a 2 C+a 3 )/ai, B = -(a 5 C+a 6 )/a 4 , A' = b 1 A+b 2 , 
B' = b 3 B + 64, and C = b 5 C, where 



Solving Eqs. (|8|)- (|T0l) . we obtain the solution of the 
Usadel equation as follows: 



C 



-030407 — 010608 + 0404040 

a 2 0407 + 010508 — 010409 



(Bl) 



ai 


= 1 - i"fB^k + - ijcf) cos 9 + (1 -+ 


- i-f B £,k + 


— ijc/, cos 9)b\ 




= h<t>(l + 6 5 )sin6»/2, 






a3 


= (1 + i~iB£k + - ijj cos 9)b 2 , 






04 


= 1 - ijB^k- + ij<f, cos 9 + (1 -+ 


■ i"fB^k- 


+ i70 cos 9)b 3 


05 


= ^(1 + 65) sin 0/2, 






06 


= (1 + «7s?fc- + H<l> cos 9)b4 , 






07 


= i70sin6>(l + 61) , 






ag 


= i70sin6>(l + 63) , 






09 


= 1 - ijB^ko + (1 + ijB^h)b 5 








= 170 sin 9(b 2 + 64) , 







8 



and 



hi 



9s~9s + »7b^+ ( 
9° S -9s- *7b£&+ 



,2ik+L 



a ik+L 



fs 



9s~ 9s ~ *7sCfc+ ' 

9°S +9s + '7Bgfc- c 2ifc-L 

9°s + 5s _ HB^k- ' 

„ik— L f+ 
e J_S 



9s + 9s~ *7sCfc- ' 
g°S + HB^kp ^ 2iknL 
9° s ~ HB(,ko 



with jf - (/o ± f s )/2. 
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